Diamagnetism in disordered graphene 
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The orbital magnetism is studied in graphene monolayer within the effective mass approximation. 
In models of short-range and long-range disorder, the magnetization is calculated with self-consistent 
Born approximation. In the zero-field limit, the susceptibility becomes highly diamagnetic around 
zero energy, while it has a long tail proportional to the inverse of the Fermi energy. We demonstrated 
how the magnetic oscillation vanishes and converges to the susceptibility, on going from a strong-field 
regime to zero-field. The behavior at zero energy is shown to be highly singular. 



I. INTRODUCTION 

The monolayer graphene has a band structure analo- 
gous to the massless relativistic particle, and its peculiar 
electronic properties have attracted much interest. Re- 
cently several experimental techniques make atomically 
thin graphene sheets available ^^i^i^ and the nature of 
this unique system is being revealed. In this paper we 
present a theoretical study on the orbital magnetism of 
graphene including the disorder effects. 

The graphene has a semi-metallic electronic structure 
where the conduction and valence bands touch at the 
Brillouin zone corners, K and K' points. Around the 
band touching point (set to e = 0), the low energy spec- 
trum has a linear dispersion analogous to the massless 
Dirac Fermion. The spectrum in a magnetic field is 
different from that in usual metals in that the Landau 
level spacing is not even but wider in lower energies, 
and is proportional to ^/B, not to B, where B is the 
magnetic field, ^ and this leads to an unusual behavior in 
the orbital magnetization. The magnetism of graphene 
was first studied as a simple model for three-dimensional 
(3D) graphite, where the susceptibility of the disorder- 
free graphene was calculated within the effective mass 
approximation. It was found that the system exhibits 
a large diamagnetism at — 0, expressed as a delta 
function of Sp at the absolute zero temperature. The 
graphene magnetism was considered again in studies on 
the graphite intercalation compounds, where the tight- 
binding model was applied for a wide range of Fermi 
energies j^i^iiSiiS. 

In the presence of the disorder, it becomes nontriv- 
ial how the magnetization behaves under this unusual 
electronic structure. Particularly, it is not clear how the 
delta-function in the susceptibility is broadened, since 
we naively suppose that the scattering is absent at e = 
where the density of states vanishes. Moreover, we do 
not know how the magnetic oscillation is destroyed by 
the disorder when we go from the high-field to the low- 
field regime, and how it converges to the zero- field limit. 

The effects of disorder on graphene under magnetic 
fields have been examined in early theoretical studies be- 
fore the experimental discovery of graphene, where the 
electronic structure,— the transport properties , 
and the de Haas- van Alphen effect^ were investigated. 



More recently the Shubnikov-de Haas oscillation was 
studied in disordered graphene,— and the spectral and 
transport properties were examined in presence of lattice 
defects under magnetic fields. 

The purpose of this paper is to calculate the mag- 
netization of disordered graphene in arbitrary magnetic 
fields, and to obtain the perspective which connects the 
high-field and zero-field limit. For the model disorder, 
we introduce the short-ranged and l ong-rang ed scatter- 
ers following the formulation in Ref. Il0l)lllfl2 . and treat 
the disorder effects within a self-consistent Born approx- 
imation (SCBA). The paper is organized as follows: In 
Sec. [H] we briefly discuss the effective mass Hamiltonian 
and the SCBA in order to make this paper self-contained 
although fully discussed previously.— The analytic dis- 
cussions of the magnetization in the zero-field limit and 
the numerical calculation for finite fields are presented in 
Sec. mil Discussions and conclusions are given in Sec. 
llYl 



II. FORMULATION 

A. Hamiltonian 

We start with the effective mass Hamiltonian in an 
ideal graphene in a magnetic field given bylS 
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Ho = T 

n 



V 



(1) 

where tt = p e A with the electron momentum operator 
p and the vector potential A — (0, Bx) in the Landau 
gauge, and 7 — \/3a7o/2 with a being the lattice constant 
and 70 the hopping integral between nearest-neighbor 
carbon atoms. A graphene is composed of a honeycomb 
network of carbon atoms, where a unit cell contains a 
pair of sublattices, denoted by A and B. The Hamil- 
tonian ID) operates on a four-components wave function 
{F^,Fj^ ,F^' ,F§'), where F^ and F§ represent the 
envelope functions at A and B sites for K point, respec- 
tively, and Ff ' and F§' for K' . 

The eigenstates are labeled by (j, n, k) with the valley 
index j = K, K', the Landau level index n = 0, ±1, . . . , 
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and the wave vector k along y directionJ^ The eigenen- 
ergy depends solely on n as e„ ~ hiOB sgn(n) i/jnf, where 
hwB = ^l/l with I = ^h/eB. The wave functions are 
written as 



^ nk 



7l 



exp(ifcy) 



/ sgn(n)(-i)0|„|_i^fe 
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exp(jfcj/) 



<P|n|,fc 
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(3) 



where 0„,fc(x) = (2"n!V^0"^^^ e-^'/2^„(z), with z = 
(x + kP)/l and -ff„ being the Hermite polynomial, and 



Cn 



sgn(n) 



1 {n = 0), 
l/x/2 (ri^O), 

(n = 0), 
(n 7^ 0). 



(4) 



For the disorder potential, we consider two simple 
models: short- and long-ranged scatterers.— The first 
is on-site potential localized at a particular A or B sites 
with a random amplitude. A scatterer on A site at Ra 
is represented as 



U{r) 



1 


zaZa 




z\z'^ 



u,S{r--RA), (5) 



and that for B site at as 



f/(r) 



where we introduced zx ~ e 




Mi(5(r-R_B), (6) 



iK 



with 



X — A and _B, and Ui = {^/?,a? /2)Ui with the on-site 
energy Ui . We assume that the scatterers are equally dis- 



tributed on A and B sites with density 



and the mean square amplitude {(uf)'^) — ((uf )^) — 

Dominant scatterers in graphenes are expected to 
have their potential range larger than the lattice con- 
stant for which inter-valley scattering is much smaller 
than intra-valley scattering. Further, realistic scatter- 
ers are likely to have the range comparable to the Fermi 
wavelengthi^i2ii2^ In the following, however, we shall as- 
sume scatterers with potential range smaller than the 
Fermi wavelength. The reason is that the results are ex- 
pected to remain qualitatively the same and further that 
actual calculations are practically possible. 

In this long-range model, a scatterer at R is expressed 
by 



C/(r) = 
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UiS{r — R). 



(7) 



We assume the scatterer density rii and the mean square 
amplitude . It was shown that the transport properties 
in the short-ranged disorder and the long-ranged one are 
qualitatively similar J^iiiii^ 



B. Self-Consistent Born Approximation (SCBA) 

We introduce the self-consistent Born approximation 
for graphene, following the formulation in Ref. fldj. The 
self-energy of the disorder-averaged Green's function 
{Ga,a') is given by 

SQ,Q'(e) ~ ^ {Ua,aiUa'^,a') {Gai,a[i£)) , (8) 

with a — [j, n, k), where (• • •) represents the average over 
the impurity configurations. 

In the short-range model, the self-energy and thus the 
averaged Green's function become diagonal with respect 
to a, and further, the self-energy is independent of a. We 
then have 



{Ga.a'{£)) = 5a,a'Ga{e), 

'-^aie) = - — ^TTITj 



(9) 
(10) 



where E(£) is the self-energy. The self-consistent equa- 
tion ([5]) is explicitly written as 



S(e) = 
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(11) 



where we introduced a cutoff function g(e) which is 1 in 
\e\ ^ £c and smoothly vanishes around e — ±ec- For 
example we can take g(e) = e"/(|e|" -I- e") with a > 
2. Further, W is the dimensionless parameter for the 
disorder strength defined as 



W 



47r7^ 



The density of states per a unit area is given by 



Pie) = 
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where — gs = 2 is the valley and spin degeneracy, 
respectively. 

In the zero-field limit, ((TT|) becomes 
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E(£) ^2W tdt 
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is - E).g(t) 



(e - S)2 - f2 ■ 
The integral is approximately written in e ^ as 



E(e) = -iy(e- I])log - 



(£-E)V' 



(14) 
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where the branch of log must be appropriately chosen. 
Then we can solve this equation analytically, 



e ~ e 



2WfL 



le 



2WTq 



(16) 



where fhiz) is the Lambert W- function, which is defined 
as the inverse function of z = ye*', and 



To = £c exp 

At £ = in particular, we have 
S(0 + iO) = 



2WJ 



-iVr, 



(17) 



(18) 



In |e| 3> Fq, S is approximately written with use of the 
expansion fL{z) ~ log(2;) — loglog(z) for \z\ 1 as, 



I](e + iO) « -2We\og 



i-KWlel 



(19) 



This can be alternatively derived from (I14p with assum- 
ing |e| ^ |S|, and thus corresponds to the Boltzmann 
limit. If W ^ 1, the states around e — Ec are completely 
mixed up with those at e = 0, as expected from the 
imaginary part of S in (|19p . To avoid this undesirable 
situation, we assume 14^ <C 1 in the following calculation. 

When the magnetic field is large enough that a Landau 
level is well separated from others, pT|) can be approxi- 
mately solved around the energy of that level. The width 
of the Landau level is estimated as 2A with 



A = V2Wnujf 



(20) 



In the long-ranged model, the self-energy and Green's 
function have off-diagonal matrix elements between 
{j,n,k) and {j,~n,k). We have 

Introducing = S"^ ± the equation becomes 

(£-S-)g(£„) 



!]+(£) 

E-(£) 
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with the same W as the short-range case Eq. ([T 
density of states per a unit area becomes 



(22) 
(23) 
The 
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Im[E+(e-HiO)-t-I]"(e-HiO)] . (24) 



In a high magnetic field such that Landau levels are 
well separated, the width of the Landau level becomes 
the same as A in (|20l) for the level 7^ 0, while it is 
\/2A for A = 0. In the weak-field hmit, S+ and 
coincide and satisfy Hi 



C. Magnetization and susceptibility 



The magnetization is defined as 



(25) 



where D,(T, fi, B) is the thermodynamic potential and /i 
is the chemical potential. By noting that the electron 
concentration A^ is given by 



A = - 



.ClfJ'/ B 

we obtain so-called Maxwell's relation, 
fdM\ _ fdN\ 



(26) 



(27) 



We write in terms of the density of states p as 

/oo 
p{e,B)f{e)de, (28) 
-00 

with /(e) = 1/(1 + e^'^"^)/''^^), and calculate M by in- 
tegrating ([27]) over /i. After a little algebra, we obtain 



M 



defie) 



de 



,d p{e'.B) 
dB ' 



(29) 



In SCBA, we evaluate this by substituting p with (flS]) 
or ([M)) depending on the type of the disorder. The mag- 
netization in a nonzero temperature is always written in 
terms of that of T = as 



M(r,M) 



de 



dm 

de 



M{0,e) 



The magnetic susceptibility is given by 

dM 
'dB 



X 



(30) 



(31) 



5=0 



taking the zero- field limit in (|29p. 



III. MAGNETIZATION IN DISORDERED 
GRAPHENES 

For the short-ranged disorder, Eqs. pT|) . (fT5|) . and 
lead to the expression 



X 
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def{e)lm 



(e-E(e)) 
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(32) 



where S is the self-energy at _B = 0. The derivation of 
this is straightforward and is presented in Appendix |^ 
In the energy range e <C e^ we can use the explicit 
form for S. At absolute zero temperature, we exe- 
cute the integral to have 



3^2 ^2 To [2WT0 



(33) 
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with 



Long-ranged, n(,= 100 



F{x) 



-Im 



(34) 



The function F(x) has the maximum at a: — with 
F(0) = 1, giving 



X(0) 



37r2 ft2 m ■ 



(35) 



In the energy range |e| ^ Fq, we use ([TOl) and obtain 



X(e-F) 



37r lepl 



(36) 



which monotonicaUy decreases as {efI increases. The be- 
havior of x{^f) can be roughly described as a long-tailed 
peak (j36p which saturates around e ~ Fq to the value 



When the disorder W becomes smaller, the peak of 
the susceptibility p3p becomes narrower and higher as 
Tq behaves as cx exp{—l/2W). The integral over e rigor- 
ously becomes —g^gss'^j'^/ (GTrfi^), as proved in Appendix 
El This is roughly verified by integrating (|36p from — £c 
to £c with the region \e\ <Tq excluded and by using Eq. 
(flT)) . Thus, in the clean limit W ^ 0, x becomes a delta 
function, 



9vg, 



6^72 



(37) 



67r ;i2 

which agrees with the result in Refs. dandS 

For the long-ranged disorder, the expression for the 
susceptibility in |e| ^ Sc becomes, 



gvgs e^7^ 



xlm- 



67r2 
1 



ft2 



rfe/(e) 
3H/ 



14- VFl0g(-£2/X2) 



(38) 



with X = e — S(e). The derivation is given in Ap- 
pendix [A) Compared with the short-ranged case ([32]) . 
we have the extra second term of the order of 0{W), 
but this gives a minor effect since W is assumed to be 
small. When 0{W'^) is neglected, the susceptibility be- 
comes just 1 — iW times as large as in the short-ranged 
disorder. Accordingly the integration of x over e weakly 
depends on W, while in ^ we again get ([37)) . 

In a strong magnetic field where the Landau levels are 
resolved, the magnetization exhibits an oscillatory behav- 
ior as a function of the Fermi energy and the magnetic 
field. The damping of the magnetic oscillation in the dis- 
order was discussed in a simple approximation where the 
scattering rate was assumed to be constant. ^'^ We calcu- 
late here the magnetization at nonzero fields in SCBA, 
since this kind of treatment is essential in investigating 
the behavior at the zero energy. We numerically eval- 
uate (|29p . in which the derivative in B is taken with a 
finite increment Ai?. Here and the following we take 
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FIG. 1: Density of states (above) and the magnetization 
at T = (below) in the long-ranged disorder with several 
strength Ws. The plot is against the Fermi energy, and the 
values are per a spin and per a valley. Vertical dashed lines 
shows the energies of the Landau level in the clean limit. 



the long-ranged disorder, and plot every quantity per a 
spin and a valley. The field amplitude B is specified by 
"■c — {^c/^bY l/^i which represents how many Lan- 
dau levels are accommodated between e = and Ec- We 
set Tic = 100 here. 

As an overview of the dependence on the disorder 
strength, we plot in Fig. [T] the density of states and the 
magnetization for several VF's at zero temperature and 
a fixed magnetic field. The density of states is basically 
equivalent to that already obtained in Ref. [Tol . but we 
present this here to demonstrate the relation to the mag- 
netization. We see that the Landau levels are separated 
more clearly in the lower energy due to the larger level 
splitting, and the magnetization exhibits an oscillation in 
the corresponding region. As W becomes larger, the os- 
cillatory part vanishes from the higher-energy side. The 
results for the short-range disorder are not shown, but 
qualitatively similar to those for the long-range disorder. 

The Landau level broadening in disordered graphene 
has also been studied for a system with lattice 
vacanciesi^ The result becomes somewhat different from 
our model in that the Landau levels around e — Q be- 
come much broader than higher levels, in accordance 
with the fact that the vacancies give rise to impurity 
states around the band touching pointj ^^i^'''^^i-^^ We do 
not have a strong scattering enhancement at £ = in the 
present effective-mass model, where the on-site energy 
of the disorder potential is assumed to be much smaller 
than the 7r-band width. 

We focus on the case of VF = 0.02 and show in Fig. 
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W=0.02 
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FIG. 2: Density of states (above) and the magnetization at 
T = (below) in the long-ranged disorder with W = 0.02, in 
several magnetic fields specified by ric = [ecI^b)^ ■ Dashed 
curves show the zero-field limit. 



[5] the plots of the density of states and oi M / B for sev- 
eral different magnetic fields. We see that the oscilla- 
tion in M terminates at a certain point and in higher e 
M/B sticks to the zero-field limit x shown as a dashed 
line. The oscillation is observable when the Landau level 
spacing tihJBWn + l — ^/n\ ~ [huj bY / {'2'\e\) is larger than 
the energy broadening at S = 0, which is 7rl4^|e| in the 
Boltzmann limit p^ . Then the condition becomes 



e > 



Huif 



V2ttW 



(39) 



In Fig. [21 the boundary is indicated by an arrow, which 
actually divides the oscillating and non-oscillating parts. 

We show in Fig. [3] the renornialized density of states 
and magnetization against e/hajg. We can see that the 
Landau level width is almost independent of the magnetic 



FIG. 3; Density of states (above) and the magnetization at 
T = (below) in the long-ranged disorder with W — 0.02, in 
several magnetic fields. The vertical and horizontal axes are 
renormalized in units of factors (xhuJB- Dashed vertical lines 
are the energies of the Landau levels in the ideal limit. 



field in this scale, as expected from (pOI in the strong-field 
limit, while each level shifts toward zero energy as B be- 
comes smaller {ric larger). From the real part of E in (|19p . 
the shift can be estimated as Ae„ ^ — M^£„ log |nc/n|. 
The amplitude of the magnetic oscillation roughly scales 
as M oc hujB oc y/B, in contrast to the behavior in 
the non-oscillating region where the relation M = xB 
is valid. This is because the gain of the total energy 
U due to the magnetic field is proportional both to the 
Landau level spacing (oc a/B) and the level degeneracy 
(oc B), which gives M dU/dB oc y/B^ The oscilla- 
tion amplitude gradually reduces as B becomes smaller, 
as the level shift causes reduction of the energy gap. 

We expect that the lowest Landau gap vanishes when 
the gap width is as small as the energy broadening at 
£ = in zero field, or 



To. 



(40) 



This is equivalent with the condition that the first Lan- 
dau level is shifted onto e = 0, or |Aei| ^ hutB, as nat- 
urally expected. To focus on this critical behavior, we 
present in Fig, d] the plots of the density of states and 
the magnetization for W — 0.1, where the condition ([^0]) 
is achieved at ric ~ 20000. We see that the structure of 
the Landau level n = still survives at ric = 100, and 
M/B deviates largely from x around this region. Here 
the magnetization at e = roughly scales as M oc ^/B. 
At ric = 10000, the gap collapses and the magnetization 
peak almost reaches that of x, and thus M (x B. If we fix 
the Fermi energy at e = and start the magnetic field 
from 0, the magnetization should exhibit the crossover 
from linear to square behavior in B. 
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FIG. 4: Density of states (above) and the magnetization at 
T — (below) in the long-ranged disorder with W = 0.1, 
in several magnetic fields. Dashed curves show the zero-field 
limit. 



between the states on a Fermi surface are proportional to 
1 /fcj?, not a constant like in our simple model. This situa- 
tion is effectively modeled in our calculation by assuming 
that the parameter W depends on £i? as W oc l/e|, in 
the long-range model. Then we expect that the suscep- 
tibility in the Boltzmann limit (|36p becomes x ^ 
From the experimental value of the mobility of monolayer 
graphene, we estimate W ^ 70 /ej^ where is measured 
in units of meV. 

The magnetization becomes highly singular at zero en- 
ergy in our model, because the energy broadening, Fq, is 
exponentially small here. In the case of charged impu- 
rities, however, the scattering rate at zero energy may 
not be small since the screening effect is strongly sup- 
pressed due to the lack of the density of states^--'— We 
need a self-consistent calculation including the screening 
and the disorder to study such a case. This is out of the 
scope of this paper and left for a future study. 

The experimental measurements of the magnetization 
of two-dimensional electron systems were performed on 
the semiconductor heterostructures, by using the super- 
conducting quantum interference device (SQUID)^^^^* or 
using the torque magnetometer j ^^i^^i^^'^^ We expect that 
the detection of the graphene magnetism is also feasible 
with those techniques. 

To summarize, we have studied the magnetization in 
graphene monolayer in presence of the disorder with the 
effective mass model and the self-consistent Born approx- 
imation. The susceptibility x{^f) has a sharp diamag- 
netic peak around zero energy even in the disorder, and a 
long tail proportional to the inverse of the Fermi energy. 
We have demonstrated that with the decrease of the mag- 
netic field, the magnetic oscillation vanishes, and M/B 
converges to x as the Landau gaps arc smeared out. 



ACKNOWLEDGMENTS 



IV. DISCUSSION AND CONCLUSION 

At nonzero temperatures, the magnetization as a func- 
tion of the Fermi energy, M{e), is smoothed in accor- 
dance with (|30p . so that fine structures smaller than fcsT 
are smeared out. This effect competes with energy broad- 
ening due to the impurity scattering, denoted here as F. 
We expect the crossover from the high-field (magnetic 
oscillation) to low-field regime (M = x^) occurs when 
either of fc^T or F exceeds the Landau level spacing Ae. 
In a usual 2D metal with a constant level spacing, it is 
known that the disorder effects can be effectively included 
as the Dingle temperature, ksTo = T/'k. The reduction 
of the magnetic oscillation in disordered graphenes was 
studied with a constant F and discussed with respect to 
the Dingle temperaturei^. 

In the realistic samples used in the experiment, dom- 
inant scatterers are supposed to be screened charged 
impurities There the scattering matrix elements 
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APPENDIX A: SUSCEPTIBILITY 

We present here the derivation of the susceptibility in 
SCBA for the short-ranged disorder Eq. (|32|) and for the 
long-ranged ([55)1 . For the short-ranged case, we obtain 
from dni) and (EH), 



9vgs 
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Let us introduce a variable X = e — E to write S as a 
function (X, B) as 



T.{e,B) = t{X,B) 



W{hujBy 



E 



2 ^ X-er. 

n— — oo 



(A2) 



Using ^ = derivative of S can be written in 

terms of those of E as 



dY.{e,B) 
dB 
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dT.{X,B) 
dX 



dt{X, B) 
'dB ■ 



(A3) 



The second-order derivative can be derived similarly as 

-1 
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Equation (|A2[) can be explicitly written as 

'HO) 



w 

E(X,i?) = —At 



(A5) 



where At = {hujBf = 2j'^eB/h, and h{x) = 2X.g(\/X) 

/{x'-t) 

When |ImX| ^ ftws, h{t) is regarded as smooth with 
respect to increment At, and we can use an approxima- 
tion 



At 



hiO) 
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h{t)dt - 



^ h{nAt) 

0<Tl<OO 

{Atf 
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h'{0) + -h'{^) 



/o 12 
where 0{At^) is neglected. Then we have 

t{X,B)-t{X,Q) = -^h'{Q){At)' 
which leads to 

at 

dB 
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With (|XT]) . and (|^ . we obtain 

xlm / de' \ 1 
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(A8) 
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where E' and X' are functions of e' . Integration in e' can 
be executed with the aid of 



de' = 
and finally obtain 



aE' 

'dX' 



dX', 



X = - 



gv9s e^7^ 
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de/(e)Im 
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X2 



(All) 
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which is ([211) • We can derive the identical equation start- 
ing from the general formula based on the linear response 
theory^Si. 

We can show that the susceptibility in the present sys- 
tem has a 'sum rule', where the integration of x{^^) over /i 
is a constant independent of the disorder strength. From 
(|A12|) . we have 



xiT,fi)d^ = 



x{0,e)de 



gvds 6272 
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By replacing the integrating variable e with X, this be- 
comes 
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gvQs e^7^ 1 
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d^^ ^ + X 
dX 



X2 
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where integration path C is a circle with an infinite radius 
with clockwise direction, and we used E ^ 0(1/X) for 
large \X\. 

The susceptibility for the long-ranged disorder ((38)) can 
be derived in a similar way to the short-ranged case, while 
the procedure is rather complicated. From ([24]) and ([29]) . 
we obtain 



X 



9v9s 



Im 



2712721^ 
1 92 , 



2 952 

We introduce a variable X^ = 
t^{X+,X-,B), with 



def{e) / de' 

) J — 00 

j:+ie',B) + ^-{e',B)) 



(A15) 



B=0 

E± and define E± = 



w{hioi 



a2 \ ^ 



X .g(e„) 



X+X- - el 

n=0 " 



n— 1 



(AI6) 

(A17) 



The derivatives of E can be written in terms of E as 



dT,' 
'dB 



A, 



aE^ 

'~dB 



(A18) 
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and 



'A'dB^~ dX^dB^^ dX^dXi dB 'dBJ' 

(A19) 

where i,j, k,l = ±, repeated indices indicate summation, 
and the matrix A is defined as 



as* 

dX^' 



(A20) 



We can calculate the derivatives of at B = in 
a similar way to the short-ranged case, and then obtain 
those for through Eqs. (|A18|) and (|A19|) . As a result, 
we have 



dB^ 



(E+ + E-) 



1 



W 



B=0 
1 



1 + a + 2/3 
2/3 



2ey 
ft 

1^2 



1-a 2 4 
where X = lim^^o = Ihirs^o AT^, and 

git) 



a ^ 2W 



(3 ^ 2W 



tdt 



tdt 



X^-t^' 

(A2 - i2)2 



1 

A3 
(A21) 



(A22) 
(A23) 



Substituting (IA2ip in (jAlSp . this becomes 



'67r2 ft2 



xlm 



1 

A^ 



1 + 3Ty 



.X(£) 

de/(£) / dX' 
- 1 /3' 



1 - a' 



(A24) 



where a' and /J' have the argument A' for A in (jA22p 
and (jA23p . and the integration in e' has been replaced 
by de' (1 + a' + f3')dX. 



In the region |e| ^ Ec, (|A22p and (|A23p can be approxi- 
mately written as a « —Wlog{—el/X^) and /3 « M^. By 
substituting them in (IA24p . we can execute the integra- 
tion in A' to obtain Eq. ([55)1 . Here the expression of 
the integrand is valid only for |e| <C Sc while the integra- 
tion runs over all e, but this is justified since the integral 
converges. 
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